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Confined polar optical phonons are studied in a semiconductor double heterostructure (SDH) by 
means of a generalization of a theory developed some years ago and based on a continuous medium 
model. The treatment considers the coupling of electro-mechanical oscillations and involves disper- 
sive phonons. This approach has provided results beyond the usually applied dielectric continuum 
models, where just the electric aspect of the oscillations is analyzed. In the previous works on 
the subject the theory included phonon dispersion within a quadratic (parabolic) approximation, 
while presently linear contributions were added by a straightforward extension of the fundamental 
equations. The generalized version of the mentioned theoretical treatment leads to a description of 
long wavelength polar optical phonons showing a closer agreement with experimental data and with 
calculations along atomistic models. This is particularly important for systems where the linear 
contribution to dispersion becomes predominant. We present a systematic derivation of the under- 
lying equations, their solutions for the bulk and SDH cases, providing us a complete description 
of the dispersive modes and the associated electron-phonon Hamiltonian. The results obtained are 
applied to the case of a EuS/PbS/EuS quantum- well. 

PACS numbers: 63.22.+m, 63.20.Dj,63.20.Kr 



I. INTRODUCTION 



Polar optical phonons can be studied along the lines of continuous media treatments for bulk semiconductors or 
semiconductor microcrystallites leading to useful results valid in the long wavelength limiti^. This kind of approach 
has been also applied in the case of semiconductor nanonostructures (Quantum- Wells, Quantum- Wires, Quantum- 
Dots, etc.) and met an acceptable degree of success in spite of the nanoscopic dimensions of the system o 4 ' 5 ' 6 i 7 i 8 i 9 i 10 . The 
so-called "dielectric continuous approach" , based upon the well-known Born-Huang equations, is frequently applied 
and allows the description of the LO, TO and IF (interface) polar optical phonons. The main limitation of this 
approach should be related to the pure electromagnetic nature of the treatment, while very often just dispersionless 
phonons are addressed. Several years ago the need for a better description of polar optical phonons was recognized, 
while theories involving dispersive phonons and discussing on equal grounds both the electric and mechanical aspects 
of the oscillations were proposed. The matching boundary conditions at the heterostructure interfaces were also 
appropriately reanalyzed, putting forward one important issue: the involved oscillations display a mixed character, 
and we cannot any longer consider the TO, LO and IF phonons as individual entities. Among the various proposed 
approaches, in Refs. pTHT^Tl3| a system of coupled differential equations was derived involving the displacement 
vector u and the electric potential if for each segment of a given heterostructure, while the boundary conditions at 
the interfaces were also deduced from the fundamental equations in a cogent way. The latter approach has been 
applied during more than a decade to different nanostructures, meeting a reasonable degree of success^i^iSiiLi^. 
The obtained phonon modes are dispersive and essentially valid near the Brillouin Zone (BZ) center, while dispersion 
is assumed quadratic (parabolic) with the wave vector k. However, it has been discussed that the sole parabolic 
dispersion is not able to describe the phonons within an acceptable approximation even if we limit ourselves to a 
limited region of the BZ (near the V point The latter issue is the main motivation for the present work. 

We develop a generalization of the equations presented in the above mentioned references allowing to introduce 
a phonon dispersion that includes both linear and quadratic terms, thus leading to a frequency dependence on k of 
the general form uj 2 = ujq + (3k + (3'k 2 for the bulk semiconductor. With this aim, the equations in the previous 
referencesiiiSii&ia are modified. Concerning the boundary conditions they need not to be altered, remaining the 
same as in the cited references and in close correspondence with the underlying equations. The new equations are 
solved for two cases: the bulk semiconductor and the semiconductor double heterostructure (SDH). The SDH is a 
well known system, in which a certain semiconductor slab, made of material "1" , is sandwiched between two slabs 
of material "2". The results are applied to an interesting case, the EuS/PbS/EuS QW, where it is evident that the 
introduced generalizations play an important role It should be emphasized that the added linear terms seem to be 
necessary for achieving reliable agreement with experiment and atomistic calculations in almost all the cases we have 
analyzed. 

The paper is organized as follows. In Sec. II we propose the new coupled differential equations involved in this 
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approach and discuss their physical meaning. Section III shows the solutions for the bulk semiconductor case. In 
Section IV the solutions for the SDH are derived, including a derivation of the three possible oscillation modes and 
their dispersion laws. Section V is devoted to normalization of the phonon states and a deduction of the electron- 
phonon interaction Hamiltonian. Finally, in Section VI a physical analysis of the EuS/PbS/EuS QW is done and 
some general conclusions are also given. 

II. FUNDAMENTAL EQUATIONS 

For the relative displacement vector u and the electric potential <p we assume the equations: 

p(u 2 - u)j.)u = aVip + pfijyV ■ u - pPffi x V x u + ip(a ■ V)u + ipb(\7 ■ u) ; (1) 

and 

VV = V ■ u with 4ira 2 = puj^(e — e^) = pt^{ur L — w\) . (2) 

In the Eqs. and J2J the same notation as in Ref. ITi1Il2TIT3I is in general applied, p is the reduced mass density 
associated to the ion couple, e (eoo) is the low frequency (high frequency) dielectric constant, cut {ojl) is the transverse 
(longitudinal) phonon frequency at the T point of the bulk semiconductor, and we assume the validity of the Lyddane- 
Sachs-Teller relation (cj£/Wj> = £o/ £ oo)- The Maxwell equations are applied in the nonretarded limit (c — > oo) and 
the quantities are considered to depend harmonically on the time (~ exp(— iujt)). Parameters 0l and (3t (with 
dimensions of velocity) are introduced to describe the quadratic (parabolic) phonon dispersion and used to fit the 
bulk semiconductor dispersion curve. They were already considered in previous works. Now we are introducing new 
parameters in the form of vectors a and b, not considered in previous works and leading to linear terms in the phonon 
dispersion law. Eqs. JTJ and (J2J) represent a system of four coupled partial differential equations assumed to describe 
the confined polar optical phonons in each segment of the semiconductor heterostructure. The solutions must be 
matched at the interfaces by applying the boundary conditions both for u and ip. These boundary conditions were 
already discussed iniLiSii^ and similar papers, and remain to be valid in the present case as can be easily proved. Let 
us recall them: 

• Continuity of u and ip at the interfaces. 

• Continuity of the normal component of the electric induction vector D — tooE + Anau. 

• Continuity of the normal component of the "stress tensor" Cy, i.e., continuity of the force flux through the 
interface. 

In our case the tensor <Jij is defined in the form: 

a H = ~P(Pl - 2/5^)(V ■ u)5ij - p(i\ {uj.i + Uij) + ip {biUj + Ojtij) . 

where the suffixes i, j = 1, 2, 3 and should not be confused with the i = \J — 1 used as a factor. Obviously, in all 
the equations above we assume an isotropic model for the semiconductor compound. The comma in Ui t j denotes the 
partial derivative with respect to xj. The four latter terms at the r.h.s. of Eq. may be written as the divergence 
of tensor cr^ and represent a certain force density which should be added to the force density —pto^u. It must be 
remarked that ov,-, strictly speaking, should not be interpreted as the usual stress tensor applied in the theory of 
elastic media. In the first place, vector u is the relative displacement between the ions (and not the displacement of 
each atom with respect to its equilibrium position). But, as has been extensively discussed in all previously published 
works on this subject, the considered terms are phenomenological in their nature and introduced with the direct 
purpose of reproducing phonon dispersion as observed in experiments and in microscopic (atomistic) calculations. In 
this connection tensor aij must be considered as a pure phenomenological quantity allowing a better description of 
phonon dispersion, and we have not attempted to relate it to the elastic constants of a continuous medium. However, 
it may be assumed to introduce in our long wavelength approach an appropriate account of what is done in atomistic 
treatments, where the interaction between and atom and its different neighbors is used for the calculations. 

We must emphasize that in some of the previously published papers the third boundary condition for the tensor 
Uij has not been applied, because it was assumed that u\g = at the interfaces. The latter boundary condition holds 
whenever the mechanical oscillations from a given segment of the heterostructure do not show significant penetration 
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into the adjacent segments. This condition can be safely used when semiconductor "1" does not vibrate mechanically 
in the range of mechanical vibrations of semiconductor "2" , and this is well fulfilled when the bandwidths of the 
optical phonons belonging to (bulk) semiconductors "1" and "2" is much smaller than the frequency gap between 
them. Microscopic (atomic) models for the optical phonons also support this approximation. In this case the second 
boundary condition reduces to the continuity of e ao dip/dN\s at each interface. 

III. BULK SEMICONDUCTOR CASE 

In the case of a bulk semiconductor the quantities u and ip should show a dependence on r in the form ~ exp(ife • r), 
and direct substitution of these solutions in Eqs. QJ and J5J leads to: 

47rm 

<p = U\ k - U ) (3) 

and 

(uj 2 - u? T + a ■ k + fS 2 T k 2 ) u = [{uj 2 l - u? T )kjk 2 - /3 2 L k + (3^k - b] (k ■ u) . (4) 
Taking u = ut + Ul where k ■ ut = 0, k ■ ul ^ 0, k x ul = 0, and k x ut ^ we are directly led to: 



(lj 2 - ui 2 L + (a + b) ■ k + l3 2 L k 2 )u L = , (5) 

and 

{uj 2 - uj^ + a ■ k + l3lk 2 )u T = . (6) 

Eq. is obtained after scalar multiplication by vector k and describes longitudinal mechanical vibrations coupled 
to the electric potential ip. Eq. ijfJJ follows after vector multiplication by k and describes transversal mechanical 
vibrations not involving an electric field. The corresponding dispersion laws are evidently implicit in the equations. 
Moreover, by applying the equation D = e(w, k)E = e^E 4- Anau, we can find the dielectric function in the form: 



e(f, fe) = u 2 - uj 2 + (q + b) ■ k + f3 2 k 2 
eoo uj 2 - w% + a ■ k + (3\k 2 ' ^ ' 

which is the expected structure of the dielectric function for this type of dispersive phonons. 

Let us now remark that, within the isotropic model we are applying, the direction of vectors a and b may be 
chosen in a suitable way. Thus, we have found convenient to take them along the vector fc, so we can finally write 
a ■ k = ak and b ■ k = bk. The sign and value of the quantities a and b should be found after a fitting of the 
bulk semiconductor dispersion law, embracing a certain percent (say, 30-40 %) of the BZ near the T point. By the 
same method parameters Pt(l) should be determined. Depending on the involved semiconductor compound we may 
be led to consider the quantity p^fi,) w ith a minus sign, and this can be easily done by the formal substitution 



IV. SEMICONDUCTOR DOUBLE HETEROSTRUCTURE CASE 



We consider the SDH grown along the z-axis, such that for — d/2 < z < d/2 we have semiconductor "1", and for 
\z\ > d/2 we have semiconductor "2". Obviously, the interfaces are located at z = ±d/2 and are formally assumed as 
infinite plane surfaces. In order to find the solutions of Eqs. JTJ and (0) we shall use the auxiliary potentials A and 
^, such that: 



u = V* + V • A with V • A = . 
Substitution of Eq. (JSJ in Eqs. and J2J), provides equations for and A having the form: 



(8) 
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V 2 * 



Pi 



-^(a + 6)-V* 



= 



(9) 



and 



V 2 A- 



u% — ui 2 A i _ A i _ _ 



= , 



(10) 



As may be seen the Eqs. ^ and IjlOfl bear a certain analogy with the corresponding equations in Ref. However, 
besides the presence of linear terms in V involving first order derivatives, we notice that Eq. 1)10(1 is now coupled to 
Eq. 10. We are thus led to equations showing some significant differences with respect to those ofii^. 

In order to find the solutions of the Eqs. 10 and (|10fl we consider the wave vector k = (k x , k v ), belonging to the 
(x, y) plane, and propose: 



ij/(r) = v(z) exp(iK • R) , A(r) = A(z) exp(zK • R) 



11) 



The exponential in Eq. Ijlljl is a consequence of the translational invariance within the (x, y) plane, this invariance 
is evidently destroyed for displacements along the z axis. We are taking R = (x, y). Vectors a and b shall be taken 
along vector n. Thus, we are led to the following equation for v(z): 



cPv 
dz 2 



Q\v = s with Q L 



oj\ — uj 2 — (a + 6)k 



(12) 



where s satisfies the equation 
form: 



- K 2 



s = 0. It is then obvious that the general solution for v(z) is given in the 



v{z) = Bexp(iQ£z) + £?'exp(— iQlz) 



Ql 



(C exp(«:z) + C ' exp(— nz)) , 



(13) 



displaying four constants as corresponds to a solution of Eq. ©. The latter solution applies to the interval — d/2 < 
z < d/2, while for \z\ > d/2 the divergent solutions should be eliminated. 

In our model of an isotropic semiconductor, we may take advantage of rotational invariance about the z axis, and, 
without lost of generality, take the y axis along vector k. Hence, we can write k = (0, k). Under the latter conditions, 
the equation for A(z) is reduced to: 



d 2 A 

dz 2 



+ Q 2 T A(z) 



ib dv 

fi\ dz 



e z = s with Q\ — 



uj^ — u) — an 



ft 



(14) 



where the components of vector s are solutions of the equation -4-^ — n 2 Si = and Si with i — 1, 2, 3 represents 

the cartesian components. Vector e z is a unit vector along the z axis. We are thus led to the following equations for 
the components of A: 



dz 2 



Q 



A: 



ib dv 

ft dz 



dz 2 



Qt 



dz 2 



Q 



A, 



The solutions of Eqs. 1)1 5H are given in the form: 
A x = B x exp(iQ T z) + B' x exp(-iQ T z) 



P 2 (Q 2 T - Ql) 



[B exp(iQ L z) - B' exp(-iQ L z)] 



(15) 



1 



[(C x - iSC) exp(Kz) + {C' x + iSC') exp(-Kz)] 



(16) 
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where S = £>/c/(/3|(Q| + k 2 )). 

A d = Bj exp(iQ T z) + B- exp(-iQ T z) + — T 



[Cj exp(nz) + Cj exp( 



-IKZ) 



with j = y, z 



(17) 



The condition V -A = leads to the equation iK,A y (z) + -j^A z (z) = and gives additional relations for the constants: 
nB y — ~QtB z , nB' y = Q T B Z , iC y — —C z , and iC y = C' z . 

Using the latter results for \P and A in Eq. (JSJ, we may obtain an explicit expression for u. We are thus led to: 



where 



X(z) 
Y{z) 

Z(z) 



u = [X(z)e x + Y(z)e y + Z(z)e z ) exp(iKy) , 
iQt [By exp(iQTz) - By exp(-iQTz)] , 

— iQr [B x exp(iQ T z) — B' x exp(— %Q T z)\ + ikT [B exp(iQ L z) + B' exp(—iQ L z)] 
[C exp(nz) + C' exp(— kz)] , 



in \B X exp(iQ T z) + B' x exp(—iQ T z)] + iQl^' [B exp(iQ L z) — B' exp(—iQ L z) 
[Cexp(Kz) — C'exp(— kz)] 



(18) 

(19) 
(20) 

(21) 



with 



bQ 



P 2 t k{Q 2 t - Ql) 



V = 1 



while the constants were redefined in a convenient way. Once we have the solutions for u, the solution for the electric 
potential is easily found by using Eq. We just write the final result: 



\B exp(iQLz) + B exp(— iQlz)] — C exp(nz) — C exp(— kz) > exp(iKy) = f{z) exp(iKy) . (22) 



[ c 



It may be seen that the X solutions represent uncoupled oscillations of a pure transversal character and not involving 
an electric potential. However, solutions for Y, Z and tp are coupled solutions. Up to this point we have not taken 
into account the system invariance under the transformation z — > —z. If the latter symmetry is used, the solutions 
can be classified according to their parity. Then, for the coupled solutions, we can write: 

Case (a): Odd Potential States (OPS) 



Y(z) 

C 
Z{z) 

c 

m 

c 



— kTSi sin(QLz) + QtS% sxd.(Qtz) + S3 sinh(Kz) , 
iQiX 'Si cos(Qlz) + 1KS2 cos(Qtz) — iSs cosIi(kz) 



47ra . 
eioo 



S lS in(Q L z) < ^ (Qi 



k(u> 



■ UJn 



-S3 sinh(Kz) 



(23) 



In the Eqs. I|23|) we have required u — at the interfaces. The solutions given are defined for — d/2 < z < d/2. For 
\z\ > d/2 we considered u = while the solution for f(z) is casted as: 



4nai 



C exp(^) 



Di sm?7 H — -5 oT-'- l 3 sinhf 



■ exp(— k\z\) 



and the continuity of the potential at the interfaces is ensured. We are using the notation: v = Kd/2, rj = Q^d/2 and 
also fi = Qrd/2. The functions Si (i — 1, 2, 3) are reported in the Appendix. The dispersion relation for this type of 
oscillations is obtained from the boundary condition ei ao df < /dz\ c i/2 = £2oodf > /dz\ ( i/2 (the notation /<(>) is useful 
to denote the function at each side of the interface) and may be written as follows: 



7 [(eioo/e2oo)^ cos 77 + vsmrj) [i/sinh^cos/i + fj, cosh v sin \j\ 
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— (/i 2 + v 2 ) [(eioo/e2oo) cosh v + sinh v] [ri/ 2 sin 77 cos // + r'/i77sin/x cos 77] , 
where the parameter 7 has the form: 

- <J^)d 2 



(24) 



7 = 



4/3» 



In the limit k — > Eq. I|24|) reduces to: sin fi cos 77 = and the zeros are simply given by: /1 = Nir and 77 
(27V + l)7r/2. 

Case (b): Even Potential States (EPS) 
In this case we have 



Y(z) 
C 
Z(z) 

c 

f(z) Aira 
C eioo 



inYS[ cos(Qlz) — iQtS^ cos(Qtz) + iS^ cosh(Kz) , 
-Q L T'S[sin{Q L z) - KS 2 'sin(Q T z) + Sg'sinh^) , 
S( cos(Q L z) - ^ < f L + cosh( K z) 



(25) 



valid for -d/2 < z < d/2. For \z\ > d/2 u = and 

S{ cos 77 



/(z) = — <7exp(i/) 
eioo 



TT U 3 



So cosh ^ 



exp(— k\z\ 



The functions S 1 / (i = 1, 2, 3) are given in the Appendix. The dispersion law, obtained in the same way than in 
the former case, reads: 



7 [(eioo/e2oo)f?sin77 — ^ cos 77] |/7 sinh cos /i — 7; cosh ^ sin 77] 



— (p 2 + v 2 ) [(eioo/e2oo) sinh v + cosh v\ [TV 2 cos 77 sin /x + Y' [ir\cos\ismr\\ (26) 

In the limit k —> Eq. I|26|) reduces to: cos [i sin 77 = and the zeros are simply given by: 77 = Nir and fi = 
(2N + 1)7t/2. Let us also remark that, both for the OPS and the EPS, the following relation holds: 



M 2 2 ( W \ (Pi 



M 2 =( + ) 1 2 +(w "+ »-M" 2 -7 . (27) 



indicating that /x and 77 are not independent quantities. 
Case (c): Uncoupled Modes 

As remarked above, these modes are related to the x component of u, and should be classified as follows. The even 
uncoupled modes (EUM) and the odd uncoupled modes (OUM), described by: 

X{z) = Ccos(Qtz) (EUM) ; X(z) = C' sin{Q T z) (OUM), (28) 
where the corresponding dispersion laws are: 

(EUM) cos M = with uj 2 = u)j* — clk — f3j*K 2 — ( N ° dd ^ T ^ ' , 

(OUM) sin M = with uj 2 = u] 2 T -aK-{? T K 2 -(^ N ™^{3 T j 2 ^ (2Q) 
where N even {N ot id ) denotes an even (odd) integer. 
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wavevector[2it / a ] 



1.0 



FIG. 1: Frequency (in units of cm *) as a function of the wave vector k (in units of 27r/ao). We display the first six modes. 

V. NORMALIZATION AND ELECTRON-PHONON HAMILTONIAN 

For the normalization of the oscillations we shall transform the classical held u into a quantum-held operator by 
means of the substitution u — > u — ua K | „, where a K , „ are second quantization operators obeying boson commutation 
rules. The classical kinetic energy for the oscillations should also be transformed into a quantum-mechanical operator 
as follows: 



where 



Wkin = t;P j u* -ud 3 r = ^pu 2 I a u ,t'r 



H, 



ph 



H p h = -^P^ 2 J (w t -u + u- vJ) d 3 r = T^pu 2 J u* -ud 



(30) 



(31) 



In Eq. (|31[1 H p h describes the free phonons hamiltonian operator, while its hermitian character is ensured by con- 
struction. Requiring the latter expression to be given in standard form for phonons of the type k , n (actually, we 
must have oj — > ui K . n ) H p h — hijj K ^ n (fij. n a K>n + one can hnally find the normalization constant C: 



^K.,n—\ T, 77 , (32) 

y bpcj K , n M K , n 

where S is a normalization area (within the (x, y) plane) and M K ,„ = J (Y* (z)Y (z) + Z* (z) Z (z))dz / \C\ 2 . The explicit 
expressions for M K: „ can be analytically derived and are different for the OPS or the EPS. They are given in the 
Appendix. 

Let us now deduce the electron-phonon interaction hamiltonian. We need to transform the electric potential into 
a quantum-mechanical operator by means of the substitution; cp — > cp — ipd K ^ n , According to Eqs. or the 
electric potential displays the general form 

A7TQ. 

ip = CF(z) exp(in ■ it) , 

where the explicit expressions for C and F(z) depend on the type of phonons (OPS or EPS). Now, the electron-phonon 
interaction hamiltonian may be derived by means of H e — P h — ^^i^pa-n ,n + ( P*^t, n)/2- Considering the Eq. ll-il'l) and 
avoiding immaterial phase factors we are led to: 
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0.0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 

wavevector [2ida a ] 



FIG. 2: Frequency (in units of cm" 1 ) as a function of the wave vector k (in units of 27r/ao). We display eleven coupled phonon 
modes, corresponding to the OPS (dotted curves) and the EPS (continuous curves) oscillations. The curves were identified by 
their behavior at k = 0. These confined phonons represent oscillations of the PbS layer for d = 3 nm. 



\ h u [F(z) K<n exp(iK-R)a K , n + h.c.] , (33) 

V D€looU) K n M K n 
k , n V ' ' 

where "h.c." stands for "hcrmitian conjugate" and the quantities M Ki „ and F K n (z) are different for the OPS and 
the EPS (and reported in the Appendix). The latter electron-phonon hamiltonian is obviously a function of R and 
z, and involves a summation over all possible phonon modes. 



VI. DISCUSSION OF THE RESULTS OBTAINED 



In order to illustrate the physical meaning of the foregoing results, we shall discuss a EuS/PbS/EuS SDH, a system 
which has received attention in the latter years. The EuS is a large-gap (E gap — 1.6 eV) magnetic semiconductor 
behaving as a ferromagnetic compound below Tq = 16.5 K, while PbS is a narrow gap (E gap — 0.3 eV) semiconductor. 
Both materials grow in the rock salt structure and show a rather good lattice matching (Aa/a ~ 0.6%). From the 
viewpoint of the electronic states the layers of PbS behave as quantum- wells and the EuS layers represent barriers^. 
These types of heterostructures are complicated for the calculation of electronic states due to non parabolicity and 
anisotropy of the band structur e) 21 ' 22 . In the current work we provide the confined polar optical modes of PbS 
for the mentioned heterostructure as a direct application of our theoretical results. Concerning the phonon states 
the europium chalcogenides have not been very well studied because of very high neutron absorption. However, 
reliable theoretical calculations for the bulk semiconductor can be found in Ref. |23[ From this reference we have 
taken the parameters = 4.9, wr(r) = 178.4 cm -1 and ojl(T) — 267.0 cm' 1 valid for EuS. For the PbS the 
corresponding parameters are: Urfi) = 65.6 cm -1 , u>l(T) = 205.6 cm -1 , ao = 5.92 x 10 -8 cm (the PbS lattice 
parameter), = 18.4, a = -2.198 x 10 18 cm/s 2 , b = -1.502 x 10 19 cm/s 2 , L = 471934 cm/s, fa = 306568V^T 
cm/s. The values of (3l(t) j °> an d b were determined by the authors after fitting the bulk PbS optical phonon dispersion 
laws 2 ^ covering approximately 30-40 % of the BZ near the T point. The values of a>T(L)(r) were also provided by the 
same reference, while comes from Ref. (see also Ref. l26l) . We can see a large difference between the UtQl) 
values for the two semiconductors. The corresponding difference for the u>l(T) frequencies is not that large, but in 
both cases (see, for instance, Fig. 4 of Ref. |23|) the band widths of bulk L or T phonons are relatively narrow. Under 



-1.0 -0.5 0.0 0.5 1.0 

z/d 



FIG. 3: Electric potential energy (in units of Wo) as a function of z. We have chosen re = 0.2 (units of 2n/ao) and ui = 191.67 
cm -1 , corresponding to an EPS. 

these circumstances, as discussed in the Section 2, the validity of the boundary conditions we have used (u\s = 0) is 
well satisfied. 

Let us first consider the so-called uncoupled modes, described by Eqs. (|28J) and l|29|l . In the Fig. 1 we can see the 
frequencies of the first six uncoupled modes as a function of the wave vector. These modes represent pure mechanical 
oscillations of transversal type and do not involve electric potentials. Their interaction with the electrons is weak and 
is determined through the deformation potential. Their physical importance should be related to the actual possibility 
of detecting these oscillations by means of optical experiments, specially by IR spectroscopy. Such information is of 
use in the investigation of the nanostructure and may be related to geometrical characteristics of the system. 

The so-called "coupled modes" are more interesting, involving a richer physical structure and are Raman active. 
These modes, as explained in the text above, are coupled to a long-range electric potential and were classified as 
OPS and EPS. In Fig. 2 we are showing the frequencies of the coupled modes (in cm -1 ) as a function of the wave 
vector k (in units of 27r/ao). We have considered a EuS/PbS/EuS SDH with d = 3 nm, and the phonons displayed in 
the figure correspond to oscillations within the PbS layer, i.e., they are polar optical phonons confined to this layer. 
The curves were obtained by numerical solution of Eq. (|24[) for the OPS and Eq. (|26|l for the EPS. In Fig. 2 eleven 
modes are displayed and denoted by LOi or TO,; according to their behavior at n = 0. The OPS modes are shown by 
dotted curves, while the EPS modes correspond to the continuous curves. Of course, curves corresponding to modes 
of different parity can cross each other. As we may see in the figure the dotted curves freely cross the continuous 
ones and, for some modes, practically show a high degree of overlapping. However, curves of the same parity do not 
cross each other and, when they become too close, anticrossing effects are seen. This is the case of the modes LO§ 
and TO5, which create an anticrossing effect near k = 0. In Fig. 2 we present the modes for a frequency interval from 
130 up to 230 cm" 1 , where the mentioned eleven modes were found. Actually, the number of possible phonon modes 
is limited by the thickness of the PbS layer, and we are probably showing more modes than those really possible. An 
important characteristic of the confined coupled phonons is the mixed nature they show. We no longer have pure TO, 
LO or IF phonons, but the dispersion curves may be predominantly TO, LO or IF for different values of the wave 
vector. 

However, at re = they recover their pure TO or LO nature, a fact that was helpful in other to label the curves. The 
regions showing a strong change in the curve slope (mainly near re = 0) represent modes where the interface character 
is more significant, involving larger electric potentials. It is interesting to remark that some curves are interrupted 
(see the modes LOi, LO2, and TOi) for large values of re. This effect is found in the numerical computations, and 
may be interpreted in terms of the applied approximation. The approximation we are using in this work has a limited 
range of validity and should provide reliable results for wave vectors just near the BZ T point. In Fig. 2 we are 
extrapolating our results well beyond this region and, therefore, we should expect our equations to fail for large wave 
vectors, especially if the frequencies are high. 
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FIG. 4: Electric potential energy (in units of Wo) as a function of z. We have chosen k = 0.2 (units of 2n/ao) and ui = 176.5 
cm -1 , corresponding to an OPS. 



In Figs. 3 and 4 we displayed the electric potential energy associated to the vibrations for fixed values of k and ui. 
We have selected two different parities: EPS for Fig. 3 and OPS for Fig. 4. These quantities were plotted in units of 

1/2 

Wq = \2 r Khe 2 (tjj\ — u)j*)d/ eiooWT S 1 ] as a function of z in the interval — d < z < d. They show an oscillatory behavior 
within the QW region and decay exponentially for values of z in the outside region. The potential energy as a function 
of z, and avoiding its y and x dependence, gives us a qualitative estimation of the electron-phonon interaction in the 
SDH. However, it should be kept in mind that the electron wave function also plays a significant role in the calculation 
of the interaction. When we study quantities like scattering rates or polaronic effects the interaction matrix elements 
must be calculated. In these cases the oscillatory character of the interaction hamiltonian could be responsible for a 
weaker interaction, but a final conclusion is possible just after effectively making the calculations. 

In conclusion, the main contribution of the present work is to provide a long wavelength treatment of confined 
polar optic al p honons for SDH, which has been improved with respect to previous treatments described in Refs. ITU 
Il2lll3lll4lll5lll6lll7yi8l in the sense it involves a better description of phonon dispersion near the BZ V point. The 
inclusion of linear terms together with the already considered quadratic ones appear to be of importance for almost 
all semiconductor compounds, leading to results of higher reliability. In the especial case of the EuS/PbS/EuS SDH 
this improvement plays an important role, and may be related to the results reported in Ref. Il9l 
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Appendix A. Some of the used functions 

In the OPS (Eqs. we used the following functions: 

Si = K cos fi sinh v -\- Qt sin [i cosh v , (Al) 
5*2 = kT sin r\ cosh v — QlY' cos 77 sinh v , (A2) 
53 = K 2 rcos^sin?/ + QtQlF' sin [i cos . (A3) 

For the EPS (Eqs. (|26p) the corresponding functions are: 
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&x = Qt cos fi sinh v — k sin p, cosh ^ , 

kT cos ?7 sinh f + Qt, T ' sin 77 cosh ^ , 
K 2 r sin jj, cos 77 + QtQ lT ' cos /i sin 77 



Si 
Si 



(A4) 
(A5) 
(A6) 



The quantity M K i n introduced in Eq. Ij32(l involves tediously long but straightforward calculations, and should be 
calculated independently for the OPS and the EPS. Here we just report final results. For the OPS we obtained: 



^ c2/t^2 2 i tV2 2\ i ^ oil 2 i 2\ , ^c2/t^'2 2 t^2 2\ S ^ n ^ 7 ? , ^ c2 / 2 2\ SU1 ^^ 

-S^F, +F V ) + -S 2 (n +v) + 1 S 1 <r V -Yu)— + -S 2 {u -H)-^- 
d $is S S 

— S\ sinh 2v + A( 2 1 2 ^ [(?7 2 r' - fi 2 T) sin 77 cos // - r)v(T' - T) sin fi cos 77] 
45iS 3 



i] 2 + i^ 2 

For the EPS we are led to: 



d(rj 2 — /i 2 ) 

[(i^ 2 r + ?7 2 r ') cosh v sin 77 + ^(r ' — T) sinh ^ cos 77] — 4S 2 S 3 sinh v cos ju . 



(A7) 



Mi E ^ = ~s 2 (rV + r'V) + \s\{f + v 2 ) ~s 2 (r'V r 2 . 2 )^ - ^s 2 V - M 2 )^ 

d $is S S 

H S 2 sinh 2z/ - - -z 1 2 , \(r] 2 T ' - p 2 T) sin fi cos 77 - 77^ ' - T) sin 77 cos u] 

2i> a(?7 2 — fi 2 ) 



4Si5 3 
?7 2 + ^ 2 



[(^ 2 T + ?7 2 r ') sinh ^ cos 77 — vq{T ' — r) cosh v sin 77] — 4S2S3 cosh v sin /i 



The other quantity is F(z), which can be casted as: 



Si sin(Q L ,z) 



exp(^) 



S 3 sinh(Kz) 



for -d/2 < x < d/2 



Sisin(77) + ^i|^l 5 3 sinh(/y) ftexp(-K\z\) for \z\ > d/2 



for the OPS. and: 



(A8) 



(A9) 



S 1 'cos(Q L z) - ^/^i^^ SJcosM^z) for -d/2 < z < d/2 

expM Is/cosfT?) - ^/ffi+ffi &'coahh/)] expf-wlzl) for Id > d/2 



(A10) 



for the EPS. 
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